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abstract 

We consider the discretization 

q{t + e) + q{t - e) - 2q{t) = sin {q{t)) , 

e > a small parameter, of the pendulum equation q" — si'a{q); in system 
form, we have the discretization 

q{t + e)-q{t) =ep(t + £), p{t + e)~p{t) = e sin {q{t)). 

of the system 

q' =p, p' = sin(q). 

The latter system of ordinary differential equations has two saddle points 
at A = (0, 0), B — (27r, 0) and near both, there exist stable and unstable man- 
ifolds. It also admits a heterocHnic orbit connecting the stationary points B 
and A parametrised by qo{t) — 4arctan (e~*) and which contains the stable 
manifold of this system at A as well as its unstable manifold at B. We prove 
that the stable manifold of the point A and the unstable manifold of the point 
B do not coincide for the discretization. More precisely, we show that the 
vertical distance between these two manifolds is exponentially small but not 
zero and in particular we give an asymptotic estimate of this distance. For 
this purpose we use a method adapted from the article of Schafke-Volkmer 
[lOj using formal series and accurate estimates of the coefficients. Our result 
is similar to that of Lazutkin et. al. [9]; our method of proof, however, is quite 
different. 

Keywords: Difference equation; Manifolds; Linear operator; Formal solu- 
tion; Gevrey asymptotic; Quasi-solution 

1 Introduction 

We consider the following difference equation 

q{t + e) + q{t - e) - 2q{t) = £^ sin {q{t)) . (1.1) 
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This second order equation is a discretization of the pendulum equation g" = 
sin((7). It is equivalent to the following system of first order difference equations 

(q{t + e)=q{t)+ep{t + e), 
\p{t + e)=p{t)+e sin {q{t)). 

which can be considered as a discretization of the system 

[p =sm(g). 

The latter system has two saddle points at A = (0, 0), S = {2tt, 0) and there exist 
stable and unstable manifolds. For the discretized equation (|1.2p and sufficiently 
small e > 0, these manifolds still exist. 

The system p.3p has [qo{t) , qQ{t)) , where qo{t) = 4arctan (e~*), as a hetero- 
clinic orbit connecting the stationary points B and A] it is a parametrisation of 
the curve p = — 2sin(g/2) and contains the stable manifold of p.3p at the point 
A as well as its unstable manifold at B. This curve, together with p = 2sin(g/2), 
separates regions with periodic orbits from regions with non-periodic orbits and 
is therefore often called a separatrix. Our purpose is study the behavior of this 
separatrix under discretization of the equation - it turns out that there is no longer 
a heteroclinic orbit for system (|1.2p and its the stable manifold at A and the un- 
stable manifold at B no longer coincide. More precisely, we want to estimate the 
distance between the stable manifold of (II. 2|) at A and the unstable manifold 
of (I1.2I1 at as a function of the parameter e. 

Lazutkin et. al. [9], Gelfreich (see also Lazutkin [7][8j) had given an asymp- 
totic estimate of the splitting angle between the manifolds. Starting from a hetero- 
clinic solution of the differential equation, they study the behavior of analytic so- 
lutions of the difference equation in the neighbourhood of its singularities t = ±§i. 

We show that the distance between these two manifolds is exponentially small 
but not zero and we give an asymptotic estimate of this distance. This result 
is similar to that of Lazutkin et. al. |9]; our method of proof, however, is quite 
different. 

We use a method adapted from the article of Schafke-Volkmer [10] using a for- 
mal power series solution and accurate estimates of the coefficients. This method 
was adapted for the logistic equation in Sellama[TT]. It turns out that the adap- 
tation of this method for the pendulum equation is more difficult than in the case 
of the logistic equation. 

We will show 

Theorem 1.1. Given any positive to, it is known that for sufficiently smal Eq > 
and all t S] — oo,to] there is exactly one one point w^^{t) = {qo{t),p^ ^(t)) on the 
stable unstable manifold having first coordinate qoit). There exist constants a ^ 0, 
such that for any positive Iq 

, ^ , s TT^- \ 47ra , , , /27rt\ _2£ /I _i£\ 
disty[wJi,[t),Wgi,) = — 2~ cosh(t) sm ^ je ^ +Ol-e = I, as e \ 0, 
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uniformly for —to <t<tQ, where dist^ (P, denotes the vertical distance of a 

point P from the unstable manifolds W~^. 

This result corresponds to the result of Lazutkin et. al. [9] as the angle between 
the manifolds at an intersection point is asymptotically equivalent to 
^T^J^ disti, , but we do not want to give any detail here. 

Our proof uses the following steps. First, we construct a formal solution for 
the difference equation (II. 1|) in the form of a power series in d = 2arsinh(e/2), 
whose coefficients are polynomials m u = tanh((it/e). This is done in section 2; 
the introduction of d is necessary because polynomials are desired as coefficients. 
Then, we give asymptotic approximations of these coefficients using appropriate 
norms on spaces of polynomials. To that purpose we introduce operators on poly- 
nomials series. In section 6 we use the truncated Laplace transform to construct 
a function which satisfies (II. ip except for an exponentially small error. The next 
and last step is to give an asymptotic estimate for the distance of some point of the 
stable manifold from the unstable manifold. A calculation shows that a = 89.0334 
and therefore 47rQ = 1118.8267 (See Remark l5.4p : the corresponding constants of 
Lazutkin have already been calculated with high precision (See Lazutkin et. al. 
[9]). A proo/ that a / as in [TO] or [TT] would be possible. Y.B. Suris [l2] had 
shown that a / 0. 



2 Formal solutions 

The purpose of this section is to find a convenient formal solution for equation 
(II. ip . First, we need some preparations. We put 

u : = tanh I ~M > 

%d{t) ■ = 4arctan | exp ( — —t 



Qd 



n=l 



for a formal solution of (jl.ll) . where d = e + X]5S=3^n^" a formal powers series 
in e to be determined. 

Remark. The linearization of equation (jl.ip at the point A gives the following 
equation 

Z{t + e) + Z{t -e)- 2Z{t) = £^Z{t). 

The parameter d is such that Z{t) = e"'^* is a solution of this equation, therefore 
£ and d are coupled by the relation d = 2arcsinh(e/2). 
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By Taylor expansion, we obtain 

+00 ^ 

qodit + e)+ qodit - e) - 2qUt) = e^"^ , (2.1) 

where - — rr^n^^HOe^"/^^" is an odd polynomial /2n-i('u) multiplied by ^/l — v?; 
(2n)! 

we find /2„_i(l) = 4/(2n)!. 

Using cos((7od) = 2n^ — 1, sin(god) = 2n\/l — u^, we can express our equation 
(II. 1|) in the form 



MT^)\j ^ 1 + Ad(r~)y ^ ^ - 2^,(tx) = /(£,n, (2.2) 

or equivalently 

^'^^^^ + ....M1'^^.L/.^ - W = /(^' ^^M) (2-3) 



cosh(d) + nsinh(d) cosh((i) — nsinh(d) 
where 



f{£,u,Aa{u)) = {2U cos [Aa{u) ^1 - u^^ + ^ sm (Ad(^) \/l - n^) ^ 

+00 

- ^/2n-l(tx)d2", 

n=l 

As M — > 1, the expressions T"*" and T~ reduce to 1, the denominators in (|2.3p 
simplify to e^'^ and hence equation (|2.3|) reduces to 

(e-'^ + e'^ - 2)^^(1) =e\2 + Ad{l)) - 4(cosh(d) - 1). 

This is equivalent to (2cosh((i) — 2 — e^)(2 + ^^(1)) = and hence we have 
necessarily e = 2sinh(ci/2) if we want a formal solution such that the coefficients 
have limits as u ^ 1. 

Theorem 2.1. (On the formal solution) If e = 2sinh(d/2), then equation 
h2.2^) has a unique formal solution of the form 

+00 

Arf(n) = ^A2„_i(n)<i2-, (2.4) 

n=l 

where A2n-i{u) are odd polynomials of degree < 2n — 1. 
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Remark: A similar formal solution was found using another method in [12j . 
Proof. We will use the Induction Principle to show that there exist unique odd 
polynomials Ai, A3, A5...A2n-i such that 



(2.5) 



k=l 



satisfy 



where 



(2.6) 



- 2Zn,d{u) - f{£,U,Zn,d{u)) 



1 - (T- 
l-n2 



(2.7) 



For n = 1, a short calculation shows that we must have Ai{u) = —\u and hence 
-\ud'^. We obtain 



Zi4{u) 



Ri{d, u) 



48 48 24^ ^ ' 



Suppose now that there exists Ai, A^, A^...A2n-i such that 

n 

Zn{d,u) = A2k-l{u)d^ 



j2k 



(2.8) 



k=l 



satisfies (12. 6p . (|2.7p . We show that there is a unique polynomial A2n+iiu) such 
that 

Zn+i{d,u) = Z„(d,n) +^2n+i(^x)d'"+' (2.9) 



satisfies (|2.6I) . We put 



ii„((i,^z) = i?2n+3(u)d'"+^ + 0((i2-+6) 



(2.10) 



where i?2n+3('w) is odd and deg(i?2n+3(w)) < 2n + 3. 

We substitute Zn+i{d,u) in equation (|2.7I1 . Using Taylor expansion, (12.9 
(|2TT0]1 and e = 2sinh(d/2), we obtain 



Zn+l,d{T 



1 - (^^ 



l-V? 

A o„.2 



n+l,d\ 



1 - (T-) 

l-n2 
,3 



2^„+i,d(n) 



- 2u" + l)A2n+i{u) + (4n^ - 4n)^2n+i(^) + 
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We notice that (|2.10p is satisfied if only if 

[(1 - u^fX^^^^iu)] ' + R2n+:i{u) = 



(2.11) 



This differential equation has a unique solution vanishing at u = without 
singularity at u = 1, namely 



A2n+l{u) 



" Jl R2n+-i{s)ds 
(1-^2)2 



dt. 



(2.12) 



We now show that this solution is an odd polynomial of u. It is clear that 
R2n+3{s)ds vanishes for i = 1 and as i?2n+3('S) is odd, it also vanishes for 
t = —1. It suffices to show that -R2n+3('S) also vanishes at t = ±1. Indeed, taking 
the limit of (|2.7I) as u — > 1 as we did for (|2.3|1 and using 



\im f(e,u,Z{d,u)) = e^Z{d,l) 



we obtain 



i?2„+3(l)d2"+^ = l^e'^ + e-'^ - 2 - £2 j z{d, 1) + 0(d2-+6). 

By our choice of e = 2sinh(d/2), we obtain i?2n+3(l)d^"+*^ = 0(d2"+6). Conse- 
quently i?2ri+3(l) = 0. As R2n+?,{u) is odd, wc also have i?2n+3(-l) = --R2n+3(1) = 
0. This proves that A2n+i{u) is an odd polynomial of degree (A2„+i(n)) < 2n + 1 
and ^2n+i(0) = 0. 

The first polynomials A2n-i{u) with n > can be calculated using Maple. 



n 


1 


2 


3 


A2n-^l{u) 




I 864 " 576 " 1 




( 319 „,5 , 185 „,3 3703 .\ 
1 2880 1152" 69120 1 





(2.13) 



Now, we ntroduce the operators C2,C,52,5 defined by 

C{Z){d,u) = \[Z{d,T+'2) + Z{d,T'l)) 

S{Z){d,u) = i(Z(d,T+i) - Z{d,T-^)) 
Ci{Z){d,u) = l{Z{d,T+) + Z{d,T-)) 
Si{Z){d,u) = \{Z{d,T+) - Z{d,T-)) 



where T+2 =r+(f,u), T'^ = T-{%u) and Z{d,u) is a formal power series in d 
whose coefhcients are polynomials. We can show that 



Ci = 252 + Id 
Si = 2SC 



(2.14) 
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and 

Ci {Q-G)= Ci (Q)Ci (G) + Si {Q)Si (G) 
SiiQ ■ G) = Ci{Q)Si{G) + Si{Q)Ci{G)) 

(2.15) 

C{Q ■ G) = C{Q)C{G) + S{Q)S{G) 

S{Q ■ G) = C{Q)S{G) + S{Q)C{G) 
if Q, G are formal power series in d whose coefficients are polynomials of u. 

3 Norms for polynomials and basis 

In this section we recall some definitions and results of [TU]. Using a certain 
suquence of polynomials, we define convenient norms on spaces of polynomials 
which satisfies some useful proprieties. We denote by 

• V the set of all polynomial whose coefRcents are complex 

• Vn the spaces of all polynomials of degree less than or equal to n 
Proposition 3.1. flO^ . We define the sequence of polynomials Tn{u) by 

To{u) = 1, Ti{u) = U, Tn+l{u) = -DTn{u) for n > 1, 

n 



where the operator D is defined by 



Then we have 



1. T+{d,u)=j:Zo^n+l{u)d\ 

2. Tn{u) has exactly degree n and hence tq{u), ...,Tn{u) form a basis ofVn, 

3. r„(tanh(z))= (;^(f)"-'(tanh(z)) 

Definition 3.2. Let p G Vn- As tq{u), ...,Tn{u) form a basis ofVn, we can write 
p £Vn as 

n 

Then we define the norms 

n 
=0 



/7r\ 

P\\n = 2^\ai\[i:) ■ (3.1) 
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Theorem 3.3. /7^ . Let n,m he positive integers and p G Vn, Q £ 'Pm- The norms 
lis. 1\) have the following properties: 

1. \\Dp\\n+l < n\\p\\n. 

2. If the constant term of p in the basis {ro,ri..,T„} is zero, we have 

IblU < \\Dp\Ui. 



3. There exists a constant M2 such that \\pq\\n+m < -^2||p||n|kiim- 



4- There is a constant M3 such that that for all n > 1, \p{u)\ < M3 (|:) 
(-l<u< 1). 



5. There is a constant M4 such that for all n > 1 with p{l) = p{—l) = 

< M4p\\n 



n-2 



4 Operators 

In this section we will use some definitions of Schafke-VolkmerfTO] and adapt their 
results on operators on polynomial series to our context. Let 

{00 
Q{d,u) = ^Q„(u)d", where Q„(n) G Vn, for aU n G N 
n=0 

By abuse of notation, let \\Q\\n = IIQnlU for a polynomial series 

00 

Q{d,u) = Y,Qn{u)d^- 

n=Q 

Definition 4.1. Let f be formal power series of z whose coefficients are complex. 
We define a linear operator f{dD) on Q by 

00 n 

f{dD)Q{d, u) = Y,{Yl fiD'Qn-i{u)) (4.1) 

n=0 i=0 

where f{z) = YX^o fi^^ "-^d Q e Q. 

By the above Definition and (1) of Proposition 13.11 we can show that 
Q{d,T+{ed,u)) = {e^p{edD)Q){d,u) for Qe Q and all G C 



8 



Thus with (|2.14p and (1) of Proposition [3TT] we obtain 

i n\n 

(4.2) 



C{Q) = cosh(f D)Q, S{Q) = sinh(f D)Q 



Ci(Q) = cosh{dD)Q, SiiQ) = smh{dD)Q 

for polynomial series Q in Q. 

Remark. According to the definition of norms in (|3.1|) . we have 

If QeQ, then dQeQ and ||(iQ|U = ^||Q||„_i for all n > 1. (4.3) 

Theorem 4.2. Lei /(z) 6e formal power series having a radius of convergence 
greater than 27r and let k be a positve integer. There is a constant K such that: If 
Q is a polynomial series having the following property 



\\Q\\n< 



for n < k 

M(n -fc)!(27r)-" for n > k 



where M is independent of n and Q £ Q then the polynomial series f{dD)Q 
satisfies 



l|/(£^^)QI|n < 



for n < k 

MK{n-k)\{2'n)"^ for n > k 



Now we define on Q the following operator 

^ = 

^ 1 
where the notation ^ means simply F{dD) with F{z) = - sinh(|). 

Lemma 4.3. For each integer k there exist a positive constant K such that: If Q 
is a polynomial series with odd Qn of degree at most n, \\Q\\n = for n < k in 
case of positive k and 

WdDQWn < M{n - k) !(27r)"" for n > max(0, k), 

where M is independent of n, then the polynomial series J^~^{Q) satisfies 

{n-k+l)\ fork<l 
\J-\Q)\\n<MK{27r)-^ { (n-l)!log(n) for k = 2 

(n - 1) ! fork>3 
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Proof. We can see easily that = vrC^^ + g{dD), where C = cosh{jdD) and 
g{z) is analytic for \z\ < 4tt, and use the proof of |10| . 

We have S = dD J = J dD, but using this relation for the inversion of S 
would give an insufficient result. Using of the formula 

2 z z 
1 = — sinh( — ) + F(z)z, where F(z) = z~'^{z — 2 sinh( — )) 

is an entire function, we obtain the relation 

Q = 2JQ + F{dD) dDQ (4.5) 
for polynomial series Q & Q. This will be essential in the proof of 

Theorem 4.4. For each integer k there exist a positive constant K such that: If 
Q is a polynomial series with odd Qn of degree at most n, \\Q\\n = for n < k in 
case of positive k, and 

||5(Q)||„ < M{n - k) !(2^)-" for n > max(0, k), 

where M is independent of n, then the polynomial series Q satisfies 



WQWn < MK{27r) 



{n-k + l)\ fork<l 
(n — 1) ! log(?i) for k = 2 
/{n-l)\ fork>3 



Proof. By the preceding theorem, we have the wanted inequalities for dDQ = 
J~^SQ in the place of Q. Here we used again UdDZHn < (n — l)||Z||„_i for any 
polynomial series Z S Q. Using theorem |4 . 21 implies the same for F{dD)dDQ with 
the entire function F of (14. 5|) As \\Z\\n < ||dZ?.Z'||„+i by theorem 13.31 we find the 
wanted inequalities (and even something better in the cases k > 2) also for JQ 
because dDJ = S. Thus formula (|4.5p yields the result □ 

In order to obtain an asymptotic approximation for the coefficients of the formal 
solution, we will need to reverse some operators. This is not possible for the 
operators S and dD on the set Q, but we can define a subset Q* of Q on which 
these operators have a right inverses. 

If we define 

Q* ■= |Q(d,n) = ^P„(n)(i", where P„(u) G P*, for all n > l| . 
where is the subspace of Vn defined by 

n 

^n — lZ]"^^^^^"'! "0 = } 

i=0 
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Then, the restrictions of the operators dD,S to Q* , denoted here by the same 
symbols 

dD : Q* ^ {l-u^)d^Q 
S : Q* ^ il-u^)d^Q 

are bijective. We denote by T the inverse of the restriction of 5 to Q* , , and we 
have 

TS = ld on Q* 
Theorem 4.5. JJ^ We consider a polynomial series 



n— 1 



r„(n)d" 



Qa{d, u) ='^an{n - 1) .y^^ 

n=2 

where an = 0{n~^) as n ^ oo with some integer k >2. Let 

4 

a := - an, 

n=l 

then the coefficients {'^{Qa)}n of T{Qa) satisfy 

{r(g,)}„ - a{n - 1) = 0((n - A;) !(2^)-") 

n 

as n —> oo for n. 

Proof. The proof of this theorem is completely analogous to that of ^0 



Theorem 4.6. J7^ Let k, l,p, q be integer with p > k and q > I. Define m as the 
minimum of k + q and l+p. then there is a constant K with the following property: 
If P and Q are polynomial series such that \\P\\n = for n < p, \\Q\\n = for 
n < q and 

ll^'lln < Mi(n- A;)!(27r)-" forn>p 
IIQIIn < M2(n-/)!(2^)-" forn>q 

then 



\PQ\\n < KMiM2{n - m) ! (27r)~" forn>p + q. 



Remark 4.7. Observe that the results of this section can also be applied, if the 
constants M are replaced by any increasing sequence {Mn)neN- In theorems 14.21 
and 14.61 the first n terms of the resulting polynomial series only depend of the first 
n terms of the given series, so the "M" in the result simply has to be replaced by 
"M„". In lemma 1^31 and theorem 14.41 the first n terms of the result depend of the 
first n + 1 given terms, so "M" in the result has to be replaced by "M„+i". 
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5 Asymptotic approximation of the coefficients of the 
formal solution 

In this section we will estimate the coefficients of the formal solution obtained 
previously (section 2). The idea is to write equation (12. 2p essentially in the form 



V{d, u)S{QiS{Q2 A)){d, u) = g{d, u, A{d, u)) , (5.1) 

where V, Qi and Q2 are known polynomials of d and u and g is a certain function 
of d,u and A involving the operators S,C and J multiplied by sufficiently high 
powers of d. 

Thanks to this equation, we will estimate the coefficients of the formal solution 
using the results of the previous section. We show that the coefficients of this 
formal solution is Gevrey-1, more precisely = 0(n ! (27r)-") . 



5.1 Rewriting of equation (12.21 1 

Consider the decomposition 



A{d,u) = U{d,u) + F{d,u) (5.2) 
where U is the initial part of A calculated before 

TT,, ^ 1 j2 /^91 3 47 \ 4 / 319 5 185 3 3703 \ g 

U{d, u) = — ud^ + u-^ u ]d^ + { H u'^ u d^ . 

^ ' ' 4 V864 576 y V 2880 1152 69120 ) 

We insert this into with (ITTill and ([2J5]) . and obtain 



2 cosh((i) • Ci (F) - 2u sinh(d) ■ Si{F) = Wq ■ F + fA d, n, F{d, u) (5.3) 
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where 



Wo = ( cosh2(d) - smh2(d)) 2 + (2u^ - \)e^ cos (u ■ \J\ - u^^ 
- 2ue^ sin (u ■ y/x - u^^ ■ \/\-u^ 

fi(d,u,F{d,u)'^ = £2 cosh^ (ci) -u2sinh2(d)) (^ ' ^^~\ sm (uWl-u'^) 

+ 2u cos (u ■ Vl - -u^) ^ cos -u) Vl - u^) + 

(^ ^^ ~ \ cos (U ■ Vl - v?^ - 2usin (|i7 • Vl-'"^)^ x 
sin u) Vl - - F{d, u) 



+00 



- (cosh2(c/) - ■u2sinh2(d)) ^/2„-i(u)c/2" - 2cosh(d)Ci(C/) 

n=l 

- 2u sinh(d)5i (C/) - 2 cosh^ (d) - ginh^ (d) j . U. 
Observ that /i has the form 



00 -. 

h{d,u,F{d,u)) = y^{d,u)+yi{d,u)J27^i^-^^rHd,uf''+ 

n=l ^ 



n=l 



(5.4) 



2n+l 



where yn{d,u),n = 1,2,3 are convergent polynomial series. 
Now, we let 

F{d,u) = Q{d,u)-G{d,u), 
J{d,u) = Qi{d,u) -SiG), 

where G is a formal power series whose the first term contains d^ and 



(5.5) 



_L ( _ 319„,6 I 1079„,4 _ J3L,,2 _l ML] rfi 
"I" ' 960 " 1728 " 2880 " 8640 ) " ' 



Qi{d,u) = {u^-l)d^ + l{l-u^)d^-^il-u^)(^lu^ + u'' + l^d^ 



(5.6) 



997 



J6L^6 , 185 4 _ Jl::^^2 j8 
2160 " ^ 432 4320 ' 
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The choice of Qi{d,u) and Q{d,u) depends in a precise way of the form of the 
equation (|5.ip and has been determined using Maple. 

Using (|5.5p . (|5.6|1 and (|2.15p . we can rewrite equation (|5.3p in the form 



WoQ + fi{d,u,Fd{u) 



Using (|2.14p . we obtain 



+ 



2cosh(d)Ci(Q) - 2u smh{d)Si{Q) Ci(G) 
2cosh(d)cSi(Q) - 2u sinh{d)Ci{Q)]Si{G). 



V-S\G) + W- SC{G) =WiG + f2{d,u, Fd{u) 
where f2(d,u, Fd{u)^ = lfi(d,u,Fd{^ 



(5.7) 



V{d, u) 
W{d, u) 

Wi{d,u) 



cosh((i)Ci(Q) - u sinh(d)5i(Q) 
cosh.{d)Si{Q) — u smh{d)Ci{Q) 

- 2cosh(d)Ci(Q) + 2u sinh(d)5i(Q) + WqQ 



The calculation of the first terms of the series Wi by Maple shows that the con- 
vergent polynomial series Wi{d,u) begins with a term containing d^^. 
Using fSSl) and ([2^5]) . we find 



5(J) = 5(Qi5(G)) = CiQi)S\G) + S{Qi)SC{G). (5.9) 



Using 



Vi{d,u) 



+ 



1 + (1 - u^)d^ + 
/1351 



2160 



193 
144 



71 

432' 
49 
60 ' 



1 



tU 



12 
11 
108 



432 J 



we obtain 



Vi-CiQi) 
Vi-S{Qi) 



QiV + W2 
QiW + Ws 



where W2{d, u) and W3{d, u) are convergent polynomials series beginning with d 
With K7\l and ([SJ]) . this implies 
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Vi • 5(Qi5(G)) = W2 ■ 52(G) + W3SC{G) + QiWiG + Q1/2 (d, n, QiG((i, u)) . 

(5.10) 

This allows us to prove the following theorem 
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Theorem 5.1. 

G{d,u) = i^— + {p+-)j—^-2-^-{pd+-)H2{d,u) 

+ 5dHi{d,u) + Sid,u) (5.11) 

where a,fi,5 are constants and the polynomial series Hq, Hi, H2, S are defined by 

00 

Ho{d,u): = Yl (n-l)!(r^)"Tn(n)d" 

n=10 ^ 

n even 

00 . 

//i(d,n): = ^(n-l)!(^)" T„(n)d" 

n=9 
n odd 

H2{d,u): = ^ri!(^)" Tn{u) d^ (5.12) 

n=9 
n odd 

and S{d, u) is a polynomial series satisfying 

\\S\\n = 0({n-3)l{27rr 



To prove this theorem we need to make some overvaluations on the coefficients 
of the polynomial series S[QiS{G)) . This will be the subject of the following 
paragraph. 

Remark: Observe that the series F, G are odd in u, even in d and beginning 
with d^. The series J is even in u, odd in d and beginning with d^^ . In the series 
F, G, A, J, the degree of the polynomial that is the coefficient of (i" is at most n — 1 
; thus the results of section 4 can still be applied and d~^F, d^^G, d~^A, d^^J £ Q 

5.2 Upper bounds for the coefficients of S(^QiS{G)) 

In this paragraph, we will use equation (I5.10p . together with the definitions of Vi 
and Qi, J and G, Wi,i = 1, 2, 3, to prove 

Lemma 5.2. 

\\^S{QiS{G))\\^ = o((n - 8) !(2^)-") as n ^ 00. 



Proof. We set 

(27r)"||cS(J)|| 

en := ; " ' , " for n > 12 (5.13) 
(n — 8) ! 
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We must show that Cn = 0{n In the sequel, we will use the following conven- 
tion: if an,n = 0, 1, .... is any sequence of positve real numbers, then 

:= max(ao, ai, ...a„) for all n > 

We have 

||5(J)||^ < e„(n-8)!(2^)-" for n > 12. (5.14) 
Using Theorem 14.41 and Remark 14.71 we obtain 

||JL <i^ie++i(n-l)!(27r)-", for n > 11. (5.15) 

where Ki denotes the constant associated with the operator S in Theorem 14.41 it 
is independent of the present context. In this proof Ki,i = 1, ..9 will always denote 
constants independent of n and the sequence e^. 
Using (|5.5p . we obtain 

||QicS(G)||,^<Kie++i(n-l)!(2^)-" for n > 11 (5.16) 



Qi{d,u 
T2{u)d' 

beginning with 1, there is a constant K2 such that 



We use (5) of Theorem 13.31 and 14.61 Since is a convergent power series 



||5(G)||„ <K2e++3(n + l)!(27r)-", for n > 9. (5.17) 

Using again Theorem 14.41 (and remark and the fact that F = G/Q where Q 
is given in (|5.6p . we obtain 

||G||„<K3e++4(^ + 2)!(2^)-" forn>8, 
||i^L <i^3e++4(^ + 2)!(2vr)-" for n > 8, 

where is a constant independent of n. 

This together with theorem 14.61 implies that there are constants K4, L such 
that for all /c > 2 

< -?^4i?/^''^(n-5)!(27r)-" forn>8A: (5.19) 

where 

J"- - S+4^n-j+4 (n-5)\ , lor n ^ lo, 

i=8 \ J ■ 



with f^t^ := foT n < 8k. 



16 



Using Theorems O] and iJ] and Wi = 0{d^°),i = 1,2,3, we obtain 

W2-S^{G)\\^< K5e+_(^{n-9)l{2TT)"^ for n > 20 (5.20) 

W3SC{G) < A6e+_e(n - 9) !(27r)-" for n > 19, (5.21) 

n 

QiWiG < Kje^_^{n - 10) !(2^)-" for n > 20, (5.22) 

n 

Qif2(d,u,Fd{u)) <ir8(l + ^-|/W+)(n-9)!(27r)-" forn > 12, (5.23) 



k>2 



Now, let us take the equation (|5.10|) 



Vi-S{J) 


< 

n 


W2-S\G)\l + 


WsSCiG) 


+ 

n 


QiWiG 


+ 

n 



Using (lOOll . (lOTTl . (lOll and (IE23|), we obtain 

llFi • 5(J)||„ < ^9(1 + e+_6 + E f /i-^) - 9)' (2vr)-'^ 

fe>2 



(5.24) 



||5(J)||„ < Kio(l + e+_e + E f /-4^) - (2^)"" (5-25) 



Since, Vi is a convergent polynomial series begins with 1 , we also have 

k>2 

Using (|5.13l) . we obtain 



"^n — 



n 



for n > 12 



(5.26) 



fc>2 



lem5.2 

Lemma 5.3. Under the condition 115. 26\) . we have Cn = 0{n~^) as n —> 00. 
Proof. Let Ki > 10 !e^2 arbitrary number. We assume that 

Ki{n+p-l)\ 



e„ < 



(n-2)!(p+ll)! 
with some p > —1, N > 16. This gives for 16 < n < 



for 12 < n < iV - 4 



(5.27) 



The first and last term of the above sum are the largest, so we can easily estimate 

n-8 

y^{i+p + 3)l{n + p-i + 3)l <{p+ ll)!(n + p-4)! . 

i=8 
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We obtain 

< Kf ^ ^-tf,r'^' , for 16 < n < iV. 
~ ^(j9 + ll)!(n-5)! - - 

In a similar way, we can prove by induction that 

- ^ (p+ 11) !(n - 5) ! - ^ (p+ 11) !(n - 2) ! " " 

Using the assumption of the lemma, we obtain 

n [p+ 11) !(n — 2) ! 

Now we choose A^o ^ 16 so large that ^"^^^^jvo^^^^^^ — then p so large 

that (|5.27|) holds for N = Nq. In a first step, our considerations imply by induction 
over that (|5.27p holds for all N and hence 

/(n + p-l)!\ 

Cr, = C — ; — as n — > oo 

V (n-2)! ; 

for this possibly large value of p. 

As Ki is arbitrary in (I5.27p . we have also shown for any p> —1 that 

/(n + p-l)!\ 
"" = ^V (n-2)! J 

implies that 

'- = ^[ (n-2)! J ■ 
Consequently the last assertion is proved for p = — 1 and we have shown 

e„ = O(n-i) 
Finally we have proved that 

||5(J)||^ = 0((n-9)!(27r)-") 

and hence that 

\\^S{J)\l = 0{{n- 8) !(27r)-") as n ^ oo 
which completes the proof of the lemma. □ 



as n ^ oo 



as n — > oo 
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5.3 Proof of theorem 15.11 

Let E := ^5(J) = S{d~^J). The polynomial series E is odd in d and its coefficients 
are odd in u. We partition it 



Tn{u) + Pn-2{n - 3) ! 



i \n-3 



2tt 



Tn-2{u) + 



(I \ n— 5 
— j Tn-A{u) + En-6{u) 



(5.28) 



for odd n > 11, where , Pn and 7„ are real numbers and also E^ have at most 
degree n for all n. For the whole series E this is equivalent to 

S{d~^J) = E = Ei + d^E2 + d'^Es + d^E 



(5.29) 



where 



E^ 



Eo 



E^i 



+00 

n=ll 

+ 00 

n=9 

+ 00 

E^ 

n=7 

+ 00 



-1) 


f- 


V27r 


-1)!( 




v2^) 









n-l 



rn(n)(i- 



n-l 



Tn(n)d" 



n-l 



n=5 



Lemma [5.21 implies that 

an = 0(n-7), = 0(n-5), 7„ = 0(n-3) and \\Ejn = 0{{n - 2) !(2^)-"). 
Applying T to (|5.29p we obtain 



^ J = T(Ei) + d^T{E2) + fi'r(£;3) + ^i'''r(s) 

To the first three summands we apply Theorem 14.51 Thus we obtain 



(5.30) 





— a{n 


-1) 




-p{n 


-1) 




-7(n 


-1) 



l( 

\27r 
V27r 
A27r 



Tn{u) 



where 



a = - On, = - /3n,7 

n=ll n=9 
n odd n odd 



0((n-7)!(2^)-") 
0((n-5)!(2^)-") 
0((n-3)!(2^)-") 

oo 

n=7 
n odd 
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To the last part of (|5.30|) we apply Theorem 14.41 and obtain 



{T{E)}n 



0((n-l)!(27r)-"log(n)), 



thus altogether 



1 



{{-/)}n - a(n-l)!^— j r„-/?(n-3)! 



2tt 



n-2 



Tn-2 



7(n — 5) ! 
Using (|4.3|) we obtain 

{J}„ - a{n-2)\ 
— 7(n — 6) ! 



27r 



n-4 



O (n-7)!(27r)-" log(n) . 



27r 
V27r 



n— 1 / 1 

r„_i-/3(n-4)!(— 



n— 3 



n— 5 



O (n-8)!(27r)-" log(n) 



Remark 5.4. T/ie asymptotic of Jn gives a good approximation of a; its suffice 
to calculate , using a formal calculation software (for example: Pari), the first 
40 terms of A{d, u) by the recurrence of Section 2 and to evaluate the highest 
coefficients of Jn to get the approximation a = 89.0334. 

Next we observe that J = QiS{G), where Qi is given in (|5.6p Using part 5. of 
Theorem 13.31 we obtain 



EA 



an I 



27r 



n+l 
T2 



+ {P-ap2iu)){n-2)l 



+(7 - PP2{u) - ap4{u)){n - 4) !^^27r 
C'f(n-6)!(2^)-"log(n ' 



27r 



n-l 



Tn-1 
T2 



T2 



(5.31) 



where 



P2{u) 



1 




1 












~4 ~ 


4 


~ ~2 




1 


n2 




1 


24 ~ 


48 


~ 432 


~ 216 



55 7 

To H Ta. 

1296 432 



Remark: Observe that the approximation (I5.3ip of the coefficients {5(G)}„ 
is polynomial. Indeed; the polynomials r„(n),n > 2 are divisible by T2{u). 

In order to find an asymptotic estimation for the coefficients of the formal 
solution, we need to apply the inverse of operator S . To this purpose, we show 
the following lemma 
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Lemma 5.5. If Hq, Hi, H2 are the polynomial series defined in I15.1S\) . Then 

• 1 If we define the operator Ci = cosh^^^^, then the polynomial series 

Ci{Ho),Ci{Hi) are converging. 

4 4 

• 2. the polynomial series S{Hq),S (Hi) are converging. 

• 3. C{Ho) = —Hq + fii{d,u), where ^i{d,u) is a convergent series. 

1 d 

• 4- •S{H2) = —HQ+^2{d,u), where H2 = d—Hi and fi2{d,u) is a convergent 

2 od 

series . 

• 5. s( — rr] = —usmh.{d) — + fis{d,u), where ^^{d,u) is a convergent 

\T2{u)J T2[U) 



series 



6. S{ — -r-T ) = ~^(^^ + l)sinh(d) — + ^Ad,u), where ^Ad,u) is a 

\T2[U)J T2{U) 

convergent series 



• 7. s( — , . ^ . . — H2 ] = T^ + U5((i, n), where UL^{d,u) is a convergent 

\smh{d)T2{u) ) - 



T2W 



series 



Proof. (1)- We have 



00 -| / 00 . \ 

rrr^ — n * V ^ — in / 



m=0 ra=10 
m even n even 



using the definition of the operator D in Proposition 13.11 we obtain 



m,n 

00 



where 



Hence 



fc=10 
k even 



1 .2vr. 



m=0 

m even 



El / 27r N rrt 1 / 27r N n 

A'^ml ^ i ^ ^ 4™m! ^ z ^ 

m,=0 m=k—8 

m even m even 
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Using 

oo 



m=0 /=0 

m even 



we find 



which implies 



00 ^ 

m=fc— 8 
m etien 



- (A;-8)!V2y ^ m!^2^ 

^ ' m=0 

m even 

^ cosh(7r/2) /-vrxfc-s 



m 



(A; -8)! V2 



This with (|5.32p imply that Ci[Hq) = fi{d,u) is convergent. For Ci{Hi), we can 

4 4 

use the same method. 

(2) - As 5 = 2S1C1, where Si = sinh ( — — ), we obtain using (1), 

44 4 \ 4 / 

S{Ho) = 2Si(Ci{Ho)) =25i(/i) 

4 \ 4 / 4 

This implies that S{Hq) is convergent. For S{Hi), we can use the same method. 

(3) - We have C{Ho) = (2C| - Id)Ho = -Hq + 2C|(Fo). This with (1) imply 

4 4 

C(ilo) = -Ho + ^ll{d,u) 

where fii{d,u) is a convergent series. 

(4) - We differentiate the equation S{Hi) = Ci with respect to d. As 

z^(^sinh(z/2)) = |cosh(z/2), 

we obtain 

S{H2) + -C{dDH,) = d-^ 
Using (3) of this lemma, we obtain 

S{H2) = ^dDHi+ fi2{d,u) = ^Ho + H2{d,u), 
where ^2{d,u) is a convergent series. 
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(5)- Using (|2.15|) and (2), (3) of this lemma we obtain 



SiHo)=S{T2iu 



C{Ho)=C[t2{u 



T2{u) 



T2{u) 



S{t2{u))c(^^^ +CiT2{u))s(^^ 

C{r2{u))c(^^^ +S{r2{u))s(^^ 
-Ho. 



This imphs 



Ho 

T2 



S{r2{u)) 



C{T2{u)Y-S{r2{n))' 



Ho + /U4(d, u) 



H 



u sinh((i) — -— + /i4 (d, n) , 

T2{U) 



where fi3{d,u), ^i[d,u) are converging series. 

(6) - The proof of (6) is similar to that of (5). 

(7) - Using (4) and (6), we obtain 



uHo 



sinh((i)T2(u) 



H, 



1 1 \ H 

-{U^ + 1) - -T2{u)]—^+ll^{d,u) 

2 2 / T2(U) 



H 







T2{u) 



+ H5id, u), 



where ii^{d,u) is a convergent series. This completes the proof of Lemma. 



Using the definition of Ho in (I5.12p and (|5.31|) . we can rewrite —S{G) in the 
form 



-S{G) = -a-^^-{P + la)^ + \aHo + X, 

d d'^T2{u) 2 T2{u) 4 



(5.33) 



where 

||X||„ = 0((n-3)!(27r)- 

Using (2) and (5) of the previous Lemma and applying also the inverse of the 
operator S in (|5.33ll . we obtain 
RE34 



^G = l(a + (/3 + ia)d^ 



1 uHo 
smh{d) T2{u) 



Ho 



1 



aH2 + 6Hi+Xi (5.34) 
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where, 5 is a constant and 



H2{d,u) = d^(d,n)= J] n!(— ) T„(n)d" (5.35) 

n=9 
n odd 



\\X,\l = 0(^(n-2)!(27r)-"j (5.36) 

, oo 

S = -V<5„, (5.37) 
vr ^ — • 

n=l 

with 5n = 0(n~^). In the expression of ^G, the term 5Hi comes from the fact that 
the series X can be written 



where 



X{d, u) = Di{d, u) + d^D2{d, u) 



D,{d,u) = ^<5„(n-l)!( T„(^x)d" 

n=8 V ^/ 

\\D2\l = 0((n-l)!(2^)- 

if we apply theorem 14.51 on the series Di{d,u) and Theorem 14.41 on Di{d,u), the 
term SHi appears in the expression of ^G. 

Since — — — - = d~^ h O(d^), we obtain 

smh(a) 6 

G{d,u) = [^ + {(3+ 3 )J ^^^'^j ' -{Pd+ -^)H2{d, u) 

+ 6dHi{d,u) +S{d,u) (5.38) 



where = o(^{n- 3)!(27r)") □ 



u H 

Observe that in d"^ — — , d^^H2 the degree of the coefficients of exceeds n. 

T2{U) 

This is due to the fact that the expressions u^^^^ — r„,+i etc., which are of degree 
n — 1, were spht. 

It is not necessary (but would not be difficult) to write down asymptotic ap- 
proximations for the coefficients of F, because equations (|5.5p and (|5.3p can be 
used. This completes the proof of the theorem 15.11 □ 

6 Functions and quasi-solutions 

So far, we have shown that equation (|2.2|) has a formal solution and we have found 
an asymptotic approximation of the coefficients of the formal solution. We will use 
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this to construct a quasi-solution, i.e. a function that satisfies equation (12. 2p except 
for some exponentially small error. To that purpose, we define the functions 



Hju) 



and 



ho{t,u) 
h2{t,u) 



j.n—1 



n=10 

n even 



(n - 1) ! 



n— 1 



2tt 
i 

2^ 



n=9 
n odd 



(n - 1) ! 



j.n—1 



(6.1) 

(6.2) 
(6.3) 
(6.4) 



n=9 
n odd 



This means that 

oo 

ho{t,u) = 



n=l 
n odd 



2^ 



'271' ^ \2tiJ 

n=0 

n even 



Using part 4. of the proposition (13. we obtain 



n=l n=l 
n odd n odd 



tanh(^) 



This is obviously the difference of two Taylor expansion and thus we can write 
ho{t,u) 

- \H2 1 + Hi— + Hq — + H^—^ 



in 



tanh (e + — ) - tanh (f - — ) 



(6.6) 



Similarly, 



hi{t,u) 



n=0 ^ 

n even 

— {Hi + Hs- + i?5 -7^7 — 

27r V ^ ''2 4 6 / 



25 



or equivalently 



— 1 r it it 

hi{t,u) = — ^tanh(e+— )+tanh(e-— ) 

i / t^ t^ t^\ 



(6.7) 



The functional equations for the trigonometric and hyperbohc functions imply 



that 



ho{t,u) 



hi{t,u) 



1 (l-n2)sin(^)cos(^) 



2- cos(X)^ + ..2^.^M 



7/2 sin^^j 



+ 



T2 



47r2 IGtt^ 



647r6 



t' 



2567r8 



t ^2 



(^i-u3)sin(^) 



+ 



47r2 167r4 



647r6 



rt' 



2567r« 



(6i 



For fixed real u the functions hk{.,u),k = 0,1,2 are analytic in |t| < p, where 
p = TT^. In the subsequent definition, we consider real values of u, < u < 1, here 
hk, k = 0,1,2 are also analytic with respect to t on the positive real axis. 
We define the functions Hkid, u), k = 0..3, by 



r+oo 



r 

Ho{d,u) : = e a ho{t,u)dt for < u < I 

Jo 

I'+oo 

Hi{d,u) : = e dhi{t,u)dt for < u < 1 

Jo 

r+oo 

H2{d,u) : = e a h2(t,u)dt for < u < 1. 

Jo 



(6.9) 



We have H2(d,u) = d——(d,u). Indeed 

od 



d-^{d,u) = J \^-e~dy-hi{t,u)dt 



e~d(hi{t,u) + t- —hi{t,u))dt 







f+OO 

= Jo ' 



e dh2{t,u)dt 



The functions TCk{d, .) are real analytic; they can be continued analytically 
to the interval — 1 < « < 1 in the following way. Choose some positive number 
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M and let Fi the path consisting of the segment from to Mi and of the ray 
t t + Mi,t > 0. Let r2 the symmetric path that could also be obtained using 
—M instead of M. Recalling (|6.6|) . we can also define 



/ e"^ tanh + — ) - / e"^ tanh - — ) 



+ MoKn), (6.10) 



where 



1 3 15 315 

Ho{d,u) := -^T2{u)d^ - ■^u{u)d'^ + ■^rfi{u)d^ - j^n{u)d^ , 

for — tanh(|:M) < n < 1, where ^ = artanh(M), because the singularities of tanh 
are i(|+n7r), n integer. As M is arbitrary, this defines the analytic continuation of 
Ti.o{d, .) for — 1 < n < 1. Similarly, the real analytic continuations of Tik, k = 1,2 
are defined. 

In the sequel, we use the operator C,S also for functions. 

Lemma 6.1. Consider the functions 7ik{d,u),k = 0..2, defined in /i6.S^) . Then, 
for — 1 < -u < 1 



1 Fork = 0, 1, 



nk{d,T^^) = -nkid,u) + fi'^{d,u), (6.11) 



where T^'i^T 2 are defined in 112. and the functions fj,^{d,u),k = 0,1, 
are analytic, beginning with d^^ , resp dP 

2. For k = 0,1, S{7ik) = fj,k{d,u) , where the functions iik{d,u),k = 0,1, 
are analytic, beginning with d^^, resp 

3. For k = 0,1, CiTLk) = —T~(-k{d,u) + Xk{d,u) , where the functions 
\k{d,u),k = 0,1, are analytic, beginning with d^^ , resp d^ . 

4. S{7i2) = —7io{d,u) + fi2{d,u), where fi2{d,u) is a analytic function, 
beginnings with d^^ . 



• 5. S\ — = — o(u + l)sinh(d) — -— + iiAid,u), where the function 
\T2{u)) 2V V 

fj.i{d,u) is analytic, beginnings with d^^ . 

Proof. (1)- For = we replace u by in ([621) • Using (fHlT)]) and C{T+^) = 
^{u) + ^d we obtain for < u < 1 

r+00 ■ 
'Hoid,T+^2)= e-ho{t,T+^)dt = —1+ + ii{d,T+^) (6.12) 
Jo 47r 
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where 

T+=/ e-3tanh(e + - + — / e"3 tanh - - + — )dt 

Jo Jo ^ 

If we substitute t + irid in the first part, t — irid in the second part we obtain 

/+oo—TTid ^ r+OO+TTid ^ 
e"dtanh(^H )dt+ / e~d tanh )dt . 

-TTid 27r' y^irf vr' 

Now, we apply Cauchy's theorem 

r+°° t / it it \ 

X+ = - / e-d(^tanh(C + — ) -tanh(e- — 



t it j"^^^ t it 

+ I e"dtanh(^H )dt - l e~dtanh(^ )dt . 

Jo ^vr Jo 2tt 

Substituting t = —ids in the second part, t = ids in the third part, we obtain 

t / it it \ 

X+ = / e-d(^tanh(e+— ) -tanh(^- — 


("^ s 
— 2id cos(s) tanh + d— 
Jo 27r 

With (102]) this imphes for < u < 1 

Ho{d, r+) = -Hoid, u) + fi+{d, u) (6.13) 

where 

n) = — / cos(s) tanh + d— )ds + r"''2 ). 

By real analytic continuation, this formula is valid for — 1 < u < 1. We use the 
same method for TLo{d^T~'2),TLi{d^T^^) and obtain for — 1 < n < 1 

no{d,T~2) = -no{d,u) + fiQ{d,u) 
ni{d,T+^2) = -ni{d,u) + fit{d,u) 
ni{d,T^^) = -ni{d,u) + fi^id,u) (6.14) 



where 



d d s 1 

liQ{d,u) = — cos(s)tanh(^- — )ds + /Uo((i,r"2) 



d ds 1 

IJ,f{d,u) = j sin(s) tanh (^ + — )ds + ;Ui((i,T+2) 



d d s 1 

fi7{d,u) = --r^ / sin(s) tanh (.^ - — )ds + /ii((i, r~2) 
4vr^ Jo 2^ 

^i(d,^z) = ^n(n)d - ^^3(^x)d' + ^^5(n)d' - Y^^7(n)d' 
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(2) - Using the definition of the operator S in (|2.13p and (1) of this Lemma, the 
result is immediate. 

(3) - The proof of (3) is similar to that of (2). 

(4) - For k = 1, we differentiate (16. lip with respect to d. Because 



9r±i 1 



dd 2 



then 



n2{d,u) ± ^(i - (r±^)2)^(d,r±5) = -n2{d,u) +dfi't{d,u) 

implies 

S{n2) = -^c((l-^x2)^)+d(/i'+(d,n)-M'f(rf,^)) 

d , 2\ dTCi , . 
= + f^2{d,u) 

= ^no{d,u) + fi2{d,u) 

where fi2{d,u) is analytic function beginings with d^°. 

(5)- Using (|2.15l) and (2), (3) of previous lemma , we obtain 

1 TYq 

= --(m^ + l)sinh(d)— -- + /U4((i, n), 

2 T2{U) 

where the function ^4,{d,u) is analytic, beginnings with d^^. This completes proof 
of the Lemma. 

In the sequel we consider mq — li 0]. 

Proposition 6.2. We have 

1. Uniformly for uq <u <1, 

oo 

no{d,u) ~ ^ (n-l)!(^)"r„(n)d" a5d\0 

n=2 

n even 

oo -11 

ni{d,u) ~ Tn{u)d'' asd\Q 

n=3 
n odd 

oo -11 

n2{d,u) ~ ^'(^j Tn{u)d'' as d\0 (6.15) 

n=3 
n odd 
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2.Fori = l..3,\^{d,u)\<Kd for uq < u < 1 {d > 0) 



Proof. The proof of this proposition is similar to that of 

With the aim of applying the results of [10], we consider 5„_2(u) = Rn{u), 
where Sn{u) is the remainder term in (|5.38|) . Then Rn{u),n is a sequence of 
polynomials of degree at most n and 

||ii„||n = o((n-5)!(27r)-ilog(n)) 

Lemma 6.3. llOj If we define 

r{t,u): = V i?„(n)- — {t £ C,\t\ < ,Uo < u < 1) 

^-^ (n — 1 ! 

n=10 ^ ^ 

r{t,u): = r(7r^u) + (t - 7r^)^(7r^u) (t > vr^ uq < u < 1) 

/•oo 

n{d,u) : = r{t,u)dt 
Jo 

then 

1. r is continuously differentiahle function on the set B of all {t,u) such that 
u satisfies < u < 1 and t is a complex number and satisfies \t\ < vr^ or 
t > vr^. The restriction of r to uq < u < 1, |t| < vr^ is twice continuously 
differentiahle. for fixed uq <u <1 the function r{t,u) is analytic in \t\ < vr^ 

2. TZ{d, u) is continuous, partially differentiahle with respect to u, has continuous 
partial derivative and 

oo 

n{d, u)r^ ^ iin(n)d" asd\Q (6.16) 

n=10 

3. \ll{d,u)\<Kd^, \^{d,u)\<Kd^ for uq<u<1 (d > 0) 

The importance of our definition of TZ lies in a certain compatibility with 
insertion of the functions T~^,T~ for u. First let 



X;ii+(n)d" = x;ii„(r+)d" 

1=10 n=ll 
oo oo 

j;i?-(n)d" = Y.^n{T-)d- 



n=10 n=ll 



We obtain a new sequences R^{u), i?„ (n) of polynomials of degree at most n. This 
follows from the relation 

oo ^ 

p{T+{d,u))=Y,^D'pi^)d' ■ (6.17) 

fc=0 
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Theorem 14.21 implies 

||i?+(u)||„ = o((n-5)!(27r)-"log(n)) 
\\Rn{u)\l = 0((n- 5) !(27r)-"log(n)) 

Therefore we can use the previous lemma for R^(u), Rn{u) and obtain functions 
7^+(^^,n), 7^~((^,^i). 

Theorem 6.4. There is a positive constant K independent of d, u such that 

|7^+((i,^i) -7^((^,^+)| < Kd^e'—, for{d>0,uo<u<l), 
\n~{d,u) -n{d,T-)\ < Kd^e'"^ for{d>0,uo<u<l). 

Proof. The proof is exactly the one of [10| . 

Definition 6.5. Let V{d,u) he a function defined for < d < do and uq < u < 1. 
We say that V[d, u) has property G if 

/•CO 

V{d,u) = / e^^ q{t,u)dt {0 < d < do, uq < u < I) 
Jo 

is the Laplace transform of some function q{t, u) that has the following properties 

1. q{t,u) is defined if Uo < u < 1 and either t is complex and \t\ < vr^ or t is 
real and t > 0, 

2. q{t,u) is analytic in \t\ < vr^ for uo < u < 1, 

3. q{t, u) restricted to < t < ir'^ or t > tt'^ is continuous and the limj^^2 q{t, u) 
exists for every uq < u < 1, 

4- there is a positve constant K such that 

\q{d,u)\ < Ke^^ for t>0,{0<d<do, uq < n < 1) 

Lemma 6.6. For uo < u < 1, we have 

L If 7ii{d,u),i = 0, 1,2 are the functions of /i6.9\) then 

d^Hk{d,u) = {l-u^)nk{d,u) + o{{l-u^)e~'^^, A; = 0, 1 

and 

d^n2{d, n) = (1 - u^)n2{d, u) + ((1 - u^)e-^^ , 
where 7ii{d,u),i = 1,2,3 have property G. 
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2. Let k be a psitive integer. IfT>i,T>2 have property G and their first terms in 
the Taylor development at d = 0, begin with d'^ then 

Vi{d,u)V2{d,u) = d''V{d,u) + o(d^e-^ 



where T>{d, u) has property G 

3. Any function 'D{d, u) analytic in a neighborhood of d = has property G if 
V{0, n) = for all u, 

4- IfTl(d,u) is defined by lemma W73\ then -—TTZ(d,u) has property G 

d'^ 

5. If'Di,'D2 have property G then so do Di + 1?2) T^i — T^2 o,nd Vi • 

6. If'D{d,u) has property G then 

\V{d,u)\<Kd {0<d<^) (6.18) 
with some constant K > independent of u. 



Proof 



1. For i = 0, 

(i)-If u > 0, we have 



POD 

d^Ho = (1 - n2) / e-Tig2{t, 
Jo 



u)dt (6.19) 



where ^ 

92{t,u) = - / / ho{s,u)dsdT 

2tt s 

g2{t,u) has a logarithmic singularity at tk{s) = {2k + l)7r^ it d i for 

(A; > 0,s > 0). it is analytic in \t\ < vr^ and lim4_,^2 g2{t,u) exists. 
If we put 

/■oo ^ 

Hoidju) = / e~dg2[t,u)dt 
Jo 

where 

\g2it,u), ift<7r2 



g2{TT^,u), ift>7r2 



then Tio{d,u) has property G and 

d^Hoid, u) = {l- u^)no{d, u) + ((1 - n2)e~^) . 
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(ii)- For —no < tx < 1, where < uq < 1„ we have 

Ho{d,u) = / e~dho{t,u)dt + 27Ti^^Res(e~dho{t,u),tk{s)j 

fc>0 

= J e~dhQ{t,u)dt + 0\^{l-u^)e~^j 

IT TT 

where 2' ^ 4 ' < u < uq, this formula coincides with the formula 

/ e dho{t,u)dt 
Jo 

and extends it by real analytic continuation for — uq < u < 0. 
This implis 

d^Ho{d,u) = (1-n^) j e-dg2{t,u)dt + 0({1 - u^)d^e--j 
we obtain 

d^Hoid^u) - {l-u^)'HQ{d,u) = {l-u^) j e-'dg2{t,u)dT 



+ 0{{l-u')e-^), 

where F is the path following the real line from infinity to (tt^, 0), then along 
the vertical line from (7r^,0) to ^tt^, tt^ tan((^)^ and finally along the line 

y = tan((^)x from ^tt^, tt^ tan((^)^ to infinity. 

Since g2{-,u) is bounded on F, then 

d'^no{d,u) - {l-u'^)no{d,u) < K{l-v?)e-^, 
where K is positive constant. Finally 

d^no{d,u) = {l-u^)nQ{d,u) +o{{l-u^)e-'i'^ for {Q<uo<u<l) 

For z = 1 we can use the same method. 
For z = 2, we use the same method with 

d^n2 = {l-u^) e-^g2{t,u)dt (6.20) 
Jo 

where ^ 

g2{t,u) = — 2V i i I h2{s,u)dsdadT. 

U — j JO Jo Jo 
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2. We assume that 'Di{d,u),'D2(d,u) have property G and their first terms in 
the Taylor development at c? = begin with d'^. Then 



poo 

Vi= e~^f{t,u)dt 
Jo 

V2= e~dg{t,u)dt 
Jo 

where f(t,u),g{t,u) are analytic in \t\ < vP' and f{t,u) = 0{t^~^),g{t,u) = 
0{t^-^). 

Vi{d,u)V2{d,u) = e~d{f ^ g){t,u)dt (6.21) 
Jo 

Since 

h{t,u) = {f*g){t,u)= [ f{t,u)g{t- s,u)ds 

Jo 

= / f{t- s,u)g{s,u)ds 
Jo 

= / fis,u)g{t-s,u)ds+ f{t- s,u)g{s,u)ds. 

Jt/2 Jt/2 

For t < TT^ , the function h{t, u) is k times differentiable with respect to t and 

h'it,u) = fit,u)g{0,u) + fiO,u)g{t,u)-f{^,u)gi^,u) 

ft 



+ / f{s,u)g'{t- s,u)ds+ f'{t- s,u)g{s,u)ds 

Jt/2 Jt/2 

= [ f{s)g'{t-s)ds+ f f'{t-s,u)g{s,u)ds 

Jt/2 Jt/2 



It/2 Jt/2 

- f{^^u)g{^,u) 

h^^\t,u) = [' f{s)g^''-'\t-s)ds+ [' f(''-'\t-s,u)g{s,u)ds 

Jt/2 Jt/2 

- J2&\l,u)g(''-'-''\l,u) 

n=0 

Observe that h^''\t,u) is contninuous on [0, 27r^[, it is analytic for \t\ < iP' . 
If we put 

^h{t,u), ift<7r2 



h{t, u) 
then 



-dU^\t,u)dt, 
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has property G and 

e~^h{t,u)dt = J e~^h{t,u)dt + 0[d''e~' 
= d^ J e^^U''\t,u)dt + o(^d''e~^j 
= d^V{d,u) + o(df'e-^^ 

where V[d,u) has property G. 

For < u < 1, the proof of (3), (4), (5) and (6) is exactly the one of [lO]. This 
proof is vahd for liQ < u < 1. 

Now we have a formal solution of the equation (|2.2p and an asymptotic estimate 
for its coefficients. With the results on the functions in the beginning of this section 
we have enough information to be able to give a precise function which satisfies 
the equation (|2.2|) with an exponentially small error as d — > 0. 

In Theorem 12.11 we found that (12. 2p has a uniquely determined formal power 
series solution 

A{d,u) = U{d,u) +Q{d,u)Gid,u) (6.22) 

where U, Q are defined in (|5.2|) , (|5.6|) and G is given by (|5.38|) . This suggest that 
we put 



Vd^ 3 / r2(n) d^ 

+ ddUi {d,u)+ d~ ^TZ{d, u) 
Aid,u) = U{d,u) +Q{d,u)g{d,u) (6.23) 

for d > 0, Mo < < 1, where TCi{d,u),i = 0..2 are defined in (|6.9|1 and TZ{d,u) is 
the function corresponding to Rn-, n = 8.10.. according to lemma [Ql Using (1) 
of proposition 16.21 and (2) of lemma 16.31 we obtain 

Q[d,u) ~ G{d,u) as d \ for every no < u < 1. (6.24) 

Consequently 

A{d,u) ~ A{d,u) as d \ for every uq < u < 1. (6.25) 

Theorem 6.7. The function 0{d,u) satisfies /i5.10\) except for an exponentially 
small error. More precisely 



Vvs(QiS{G)^-W2-s\g)-W3SC{g)-QiWig-Qif2(d,u,G^ 



< Kd^(l-u^)e-^ 



uniformly for {uq < u < 1,0 < d < do,) , where K is a constant independent of d 
and u. 
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In the proof of this Theorem the functions Vi{d,u),i = 1.2.. have property G. 
Proof. We set 

j'{d,u) = Vvs(^QiS{g)^-W2-s^{g)-w^sc{g)-QiWig-Qif2{d,u,gy (6.26) 

Using (2), (4), (5) of Lemma 16.11 and (I6.23|) . we obtain 

S{g) = -(ai sinh(d)— - sinh((i) - ^)') Wo + (1 - u'^)d^Vi{d, u) + d-'^S{n) 

\ T2 2 2 / 

where Vi{d,u) has property G and 

a , ^ a, 

= ^+Pd 
d 

This with lemma [HTT] imply 

Vis(QiS{g)^ = d^V2{d,u)no + {i-u^)d^^v^{d,u) + d~^ViS{Qi)cs{n) 
+ d-^ViC{Qi)s^{n), 

where V2{d,u),V^[d,u) have property G. Because 

C{Qi) = {l-u^)d^{l + 0{d^)), 
S{Qi) = {l-u^)d^{u + 0{d^)), 
Vi = l + 0{d^) 

it is sufficient to apply Theorem 16.41 (4) of Lemma WM for IZ and (1) of Lemma 
16.61 for TLq we obtain 

ViS{QiS{g)^ = d^V2{d,u)rLQ + {l-u^)d^Vi{d,u) + o{{l-u^)d^e-^^ 

where Vi{d,u) has property G. With (1) of lemma 16.61 this implies 

VxS(QiS{g)) = {l-u^)d^V5{d,u) + o{{l-u^)d^e-^^ (6.27) 

Using the same method for the terms W2 • S'^{g)^ W'iSC{g) and QiWig, we 
obtain 

W2-S\g) = {l-u^)d^^Ve{d,u) + o(^{l-u^)d^'^e-^y 
WsSCig) = {l-u^)d^V7{d,u) + o(^{l-u^)d^^e-^y 
QiWig = {l-u^)d^Vs{d,u) + o(^{l-u^)d^^e-^y (6.28) 
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To study the term Qif2(d,u,Q] , we first treat and . 



2 

g^id, u) = al%nl + Plnl + 5'^d^nl + d-*n^ - 2ai/3i— ^0^2 + 2ai5— ?^oWl 

7"! T2 T2 

+ 2d-^ai—non-2Pi6dnin2-2d-'^PiTzn2 + 2d-^6nni (6.29) 

Using (1), (2) and (4) of the lemme WM we obtain 

g{d, uf = d'^Vg{d, u) + (d^e^^^ ) (6.30) 
With (2) of lemma 16.61 this implies 

g{d, uf = g{d, u)g{d, uf = d^Vio{d, u) + (d^e-^^ . (6.31) 
We can rewrite 

f2{d, u, g) = void, u) + yi{d, u)g^h^i {d, u, g^) + y2{d, u)d^g'^fi^2{d, u, g^) (6.32) 
where /i,2 and yj, i = 0, 1, 2 are analytic . 
Lemma 6.8. For i = 1,2, 

e-^fi{t,u)dt + o(^d^e-^j for{d > 0,uo < u < I) (6.33) 

where fi{.,u) is analytic in \t\ < vr^ and continuous in [0, 7r^[ and [7r^,oo[ 
Proof. Using (|6.30l) . we can rewrite 



g'^{d,u)=l e~^g{t,u)dt + 0{d^e''^ 



where g{.,u) is analytic in |t| < tt^, it is continuous in [0,7r^[ and [7r^,oo[, twice 
differentiable in \t\ < vr^. We obtain 



fi,i{d,u,g^) = ^fn,i{d^u)g'^ 



-.2n 

n=l 

where 

/„,i((i,n)= / e~dLpn.i{t,u)dt 
Jo 

and ifn,i{t,u) are entire functions. 

Using the proof of theorem 5.1 from we find 



e~-dMt,u)dt + 0[d^e-'r'j (6.34) 
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where the series 

oo 

* g*''){t, ufj , g*"" = g* ... * g,n times 

n=l 

is uniformly convergent to a function fi{t,u) analytic in \t\ < vr^ and continuous 
in [0, 7r^[ and [7r^,oo[, satisfies also 

\fi{t,u)\ < Kexp{Kt) for t>0,uo <u< 1. 

Then 

fi4d,u,g^) =yi{d,u) + o(d^e-^^ foi{d >0,uo<u< 1), 

where yi{d,u),i = 1,2 have property G. 

This lemma with (|6.32|) and (2) of lemma 16.61 imply 

Q1/2 (d, u, e?) = d\l - u^)Vio{d, u) + (^d^l - u^)e-^^ (6.35) 
Combining (lOTll . (lOHIl and (IOBI) . we find 

T{d,u) = d^{l-u'^)V{d,u)+lC{d,u) 
where V{d, u) has property G and 

\lC{d,u)\<Kd^{l-u^)de'^ for (0 < d < 4, uq < u < 1). (6.36) 

(■00 J 

T{d, n) = (1 - n^)^^ / e-dq{t, u)dt + /C(d, u) 
Jo 

where q{t,u) is analytic in |t| < tt^, it is continuous on [0,7r^[ and ]7r^,oo[, has a 
limit as t — > TT^ for every uq < u < 1 and satisfies 

\q{t,u)\ < Ke^\ (6.37) 

with a constant K independent of u. If q{t,u) = X]^o9"(^)^" power series 

of q{t, u) near f = 0, Watson's lemma with (I6.36P imply 

00 

J^{d, u) ~ n!(l — v?)qn{u)d^~^'^ as d \ for every uq < u < 1 . 

n=0 

On the other hand because of its definition 

?2 



Hence 



T{d,u) ~ Vi-s[QiS{g)'j -W2-s\G)-w^sc{g) 
- QiWig-Qif2{d,u,g] =Q + {)d + --- 
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since the formal series G satisfies (|5.1Q|1 . This means that all g„ = 0. Thus we 
obtain for uq < ''^ < 1 with (I6.37P 



/»00 /'OO 

{l-u^)(f I e~^\q{t,u)\dt < {l-u^)d^ e-^Ke^Ut 

Jo Jn^ 

< K{l-u^)d^e-'^ (0<d<do) 

and thus 



,2 



\T{d, u) \<K{1- u^)d-^e'^ {0 < d < do) . 
We have proved that 

Vi-s(QiS{g)'^ - W2-s\g)-W3SC{g)-QiWig -Qif2[d,u,g') 

< K{l-u^)d^e-'^, 

ioi < d < do, uo < u < 1, i.e. Q{d, u) is a quasi-solution of (|5.10|) on this interval. 
This implies that the function A{d,u) defined in (|6.23l) is a quasi-solution of (12. 2p . 
more precisely 



\ ^^y ^(d,r+) + a/ \ / A{d, T- ) - 2A{d, u)-f {e, A{d, u)) 
1 — V 1 — tt 



,2 



< Kde'^, (6.3^ 



7 Distance Between Points of Manifolds 

Clearly, if Qeit) is an exact solution of the difference equation (ll.ljl . then {qe{t),Pe{t)), 
where Pe{t) = ^ {qe{t) — qeit — e)) , is an exact solution of the the system (|1.2p . In 
the introduction, we have mentioned that the stable manifold W~ of this system 
at ^ = (0, 0) is parametrized by i — > {q~ (t) , p~ {t)) and the unstable manifold 
of p.2p at i? = (27r, 0) is parametrized by t ^ {qf{t),pf{t)), where {q~ (t) , pj (t)) 
is an exact solution of (|1.2p and qf{t) = 27r — q~{—t),pf{t) = pj{—t + e). 

In the previous section, we have constructed a quasi-solution A{d, u) for equa- 
tion (12.21) . i.e. it satisfies this equation with an exponentially small error. We 
denote by Wg, Wu the manifolds close to W~ respectively parametrized by 1 1— > 
{^-{t),ip-{t)) respectively t </?+(*)), where = ^/l - u(tfA{d,u{t)) + 

qodit) and U{t) = 27r - ^-(-t), ip-{t) = ^(^-(0 - ^-(i " e)) , if+it) = ^ {C+{t) - 

— e)). Here and in the sequel, we often omit to indicate the dependence with 
respect to e for the sake of simplicity of notation. 

We will first show that the vertical distance between some point of the 

stable manifold W~ and the manifold Wg is exponentially small. For this purpose. 
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we consider the sequence Z„ = {qn,Pn) on the stable manifold Wg , defined by 

Zn,. = ('-^'] = ^f^") = ( ^" + - = 1'2,3,.... (7.1) 

VPn+l/ \PnJ \Pn+ £Sm.[qn)j 

There is a sequence t„ such that 

e_(t„) = qn = ^l^Mtnf A{d,u{tn)) + qod{tn)- (7.2) 

The vertical projection of the point {qn,Pn) on the manifold Wg is the point 
{qn,bn), where 6„ = v?_(t„) = g{qn) := (/^-(Cn^n)) for n = 1,2,3, ... . We denote 
by A£(n) the vertical distance between the point {qn,Pn) and the manifold Wg. 
Then 

A,{n)=pn-bn. (7.3) 
For the i?i)-image of the point {qn, bn) on Ws, we find 

K+i) \bn) \bn + esui{qn 

i-{tn) +ebn+l 

e-i(e-(tn) -e-(in-e)) +esin (^-(tn)) 

2e-(in) - ^-(tn - e) + sin (^-(tn)) 
e-i(e-(t„) -e-(in-e)) +esin (^-(tn)) 



(7.4) 



As the function ^_ satisfies equation (II. ip except for an exponentially small error 
because of (|6.38|) and (|7.2|) . we have 

= [tn + e)+ e„+i (e) 

= y?_(t„ + e) + -e„+i(e), (7.5) 

where |e„4.i(e)| < ii'eexp( ) with some K independent of e and n. Since 

g {i- {tn + e)^ = ^- {tn + e) , we have 

bn+i = g(^i-{tn + e)^ +^en+i{e) 

= g[qn+i- e.n+i{e)^ +^en+i{e) (7.6) 

On the other hand, using (17.31) and the definition of (j), we obtain 

qn+i\ _ qn \ _ fqn + epn + e sm{qn) - e/^£{n) 

bn+ij \Pn - Ae(n)y V Pn + esin(g„) - As{n) 

^n+i-eA (n)Y 

Pn+i-A,{n)J ^ ^ 
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With (I73|) and this imphes 

Ae(n + 1) = Pn+l - bn+l = Pn+l - g{qn+l) 

= As{n) +bn+i- g(qn+i+eAs{n] 



Using Taylor expansion we obtain 

A,(n + 1) = {l-eg'{en+i))Ae{n) + ^{l-eg'{en+i))en+i{e) (7.8) 

where - e„+i(e) < 9n+i < Qn+i + eAe(n). 

Now g is e-close to the curve p = —2sm{q/2), hence 

ff'(0„+i) = -cos(^)+0(d). 

Thus given any positive fj, < ir, there is a positive constant c such that for all 
qi < all n and sufficiently small d, 

l-eg'{en+i) >l + ec . 

It is now convenient to write (|7.8p in the form 

A,(n) = (1 - eg'{en+i)y^As{n + 1) - ^e„+i(e) 

As A^{n) — > as n — > oo, this implies that that there is a positive constant K 
such that 



\Ae{n)\ < Ke-V^(l + ec) 



k=0 



Consequently 



In particular 



A,(n) = O0exp(-^)). (7.9) 



dist,((gi,pi),VF,) = O0exp(-^)^, (7.10) 

where (<?i,pi) is any point on the stable manifold W~ , provided qi < fJ- < it; here 
dist„ denotes the vertical distance. 

The estimate (|7.10|) can be extended to any fj, < 2tt and a starting point {qi,pi) 
with (/i < /U in the following way. The relation (|7.8p remains valid, only now we 
just have the existence of some constant c > such that 1 — eg'{9n+i) > 1 — ec 
for all n. As system (|1.2p can be regarded as a one-step numerical method for the 
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system (|1.3|) of differential equations and the starting point is at a distance 0{e) 
of its solution (go(i), 9o(i)), results on the convergence of one-step methods can be 
applied and yield that qk = qoiti + {k — l)e) + 0{e), where ^0(^1) = Qi, uniformly 
for integer 1 < k < L/e, where L is any positive constant. We choose L such 
qo{ti + L) < /i/2 < vr. Repeated application of (|7.8p now gives 



I A,(n) -A,{n+ [L/e])\ < Ke'^ ^ (1 - ec)-'= < 



-e £ 



ce 

k=0 



To the quantity Ae(n-|- [-^^/e]), inequality (|7.9p can be applied, because qn+[L/e] ^ 
/i/2 < vr. Thus (|7.9p and hence also (|7.10ll remain valid also uniformly for < 
Qi < Ml provided /i < In. 

The analogous reasoning applies to the vertical distance of a point {qi,pi) on 
the unstable manifold from the manifold Wu and yields 

dist^iiqi, Pi),Wu) = oQexp ( - (7.11) 

Another method to obtain (|7.11|) consists in using (|7.10l) and symmetry. 

Now we will estimate the vertical distance between the two manifolds Wg and 
Wu- As the quasi-solution A{d,u) is defined for —l<u=: tanh(t) < 1, we can 
define 

A'^{d,u) = -A{d,-u), 



C+{t) = ^l-u{tyA+{d,uit))+27T-qod{-t)=27T-^.i-t),{7.12) 

D,{t) = ^+{t) - C^{t) for-^<t<^ . (7.13) 
Using (|7.12p and the definition of ^_(t) we find 

D,{t) = -^l-u{tf[A{d, u{t)) + A{d, -u{t))^ - qod{t) - qod{-t) + 2tt 

= - Y^l - u{t f (^A{d, u{t)) + A{d, -u{t))'^ (7.14) 
With dSS]) and (lOal) this implies 



Ds{t) = -\Jl- u{tfQ{d, u) ai^ (no{d, u) - Hoid, -u)j - 
Pi {H2{d, u) + n2id, -n)) + 5d(niid, u) + Hiid, -u) 

where Q{d, u) is defined in (|5.6|1 . ai, /3i are defind in ()6.27|1 and Ti.i{d, u) are defined 
in (16. 9p and can be continued analytically to — 1 < < 1 as in (I6.10p . 
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Using the fact that the functions uho{s,u), hi{s,u), h2{s,u) in (|6.9p are odd, 
we can apply the residue theorem and obtain for — 1 < u < 1 that 



Hi{d,u) +Hi{d,-u) 



27riRes^e <ihi{s,u)^Sk(t) 

lm{sk{t))<0 

— 27ri Res^e~^/ij(s, u), Sfc(t)^ , i = 0,1,2, 

/m(sfc(t))>0 



where Sfc(t) = vr^ ± ^ i + Ikn'^ for A: > 0. We obtain 
Res(^e~^ho{s,u),Skit) 
Res( e^d/ii(s, u), Sfc(t) 



Res(e d/i2(s, u), Sfc(t) 
and hence 



(/c + l)7r2 27rtz 
1 1 =F 

{k + l)7r2 27rtz 
±— e dee 

{k + l)7r2 27rti 

i ~ T 

±-—(TTei^2dt)e dee 



(7.15) 



where 



(pi{t,e) = 2'Ka 
+ 



sinh 



dt 



TTa 



cosh ( f 



cos 



+ t/cosh ^ 

(?) 



sm 



(?) 



As a consequence of (|7.15p and (|7.12|) . we obtain immediately that 

^±(0) = 7r + 0(e-2e-"'/^) . 



(7.16) 



(7.17) 



Now, let us take a point (^+(t), (/?+(t)) on the manifold Wu- We suppose that 
the point (^_(fi), (/9_(ti)) is its vertical projection on the manifold Ws. We will 
evaluate the vertical distance between these two points 

dist,(t) = if+it) - ^_(ti) = if+it) - g{C+{t)), (7.18) 

where g{x) = (^^Z^(x)) . Thus by (ITTSll 

disK{t) = ip+it)-g{C-it) + Deit)). (7.19) 
Using Taylor expansion, we find 

dist^it) = ip+{t) - g{C-{t)) - D,{t)g' {v{t)), (7.20) 
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where < r]{t) < C-{t) + De{t). Here 

Vii) = Qodit) + 0{d), hence 



cos 



' —r- 



V 2 



j+0(d) 



tanh (^y^ +0((i). 



(7.21) 
(7.22) 



As g{£.-{t)) = 95-(i), this yields 



dist^(t) = ^+{t)-ip.{t)-Deit)gL{r]{t)), 

= ^(e+(t) - C+(t - e)) - -^iUt) -Ut- e)) - D,{t)g'_{v{t)), 
1 



De{t)-D,{t-e))-De{t)g'{r]it)). 



Now formula (|7.15|) applies and we obtain 

De{t-e) 



Deit) + ^(l)2it,e)ew(-^] +0{e~^ 



where 



(7.23) 



(7.24) 



'/>2(i,e) 



-27ra 



cosh 



dt 



+ 1 



tanh 



dt 



cosh ( Y 



cosh ( Y 



sm 



+ vra 



tanh ( f 



Lcosh(f) 
With (17:23]) and !n\22^ this imphes 



cos 



dist^,(t) = 



62(t,e) +tanh (^^) (/>i(t, e) 



e d + 0( -e d 
a 



4 4 

Consequently, for — - < t < - 

47rQ; 

dist,;(t) = — 7T- cosh(t) sin 



as e \, 0. 



(7.25) 



(7.26) 



(7.27) 



Combining this result with (|7.10p and (|7.11l) . we finally obtain our main result 
theorem 1.1, because ^±(t) = qo{t) + C'(e^) uniformly with respect to t on any 
finite interval. 
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